Constraints on relaxation rates for iV-level quantum systems 
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We study the constraints imposed on the population and phase relaxation rates by the physical 
requirement of completely positive evolution for open iV-level systems. The Lindblad operators that 
govern the evolution of the system are expressed in terms of observable relaxation rates, explicit 
formulas for the decoherence rates due to population relaxation are derived, and it is shown that 
there are additional, non-trivial constraints on the pure dephasing rates for N > 2. Explicit, 
experimentally testable inequality constraints for the decoherence rates are derived for three and 
four-level systems, and the implications of the results are discussed for generic ladder, A and V 
systems and transitions between degenerate energy levels. 
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I. INTRODUCTION 

Understanding the dynamics of open systems is crucial 
in many areas of physics including quantum optics Q, , 
quantum measurement theory |3(, quantum state diffu- 
sion Q, quantum chaos p|, quantum information pro- 
cessing and quantum control 0, H, T 
many efforts to shed light on these issues 
many important questions remain. 

For instance, it was recognized early by Kraus |13| . 
Lindblad [LH fl5| . and Gorini, Kossakowki and Sudar- 
shan [l6| that the dynamical evolution of an open system 
must be completely positive [24| to ensure that the state 
of the open system remains physically valid at all times. 
Unfortunately, if relaxation rates are introduced ad-hoc 
based on a phenomenological description of the system, 
the resulting equations often do not satisfy this condition. 
For example, the Agarwal/Redfield equations of motion 
for a damped harmonic oscillator have been shown to 
violate complete and even simple positivity for certain 
initial conditions |17|. Although such master equations 
may provide physical solutions in some cases, serious in- 
consistencies such as negative or imaginary probabilities, 
unbounded solutions and other problems may arise. 

For two-level systems the implications of the complete 
positivity requirement have been studied extensively in 
the literature, for example by Gorini et al. who first 
showed that there are constraints on the relaxation rates 
in the weak coupling limit 0, 0] , Lendi who provided a 
comprehensive in-depth analysis of the dissipative opti- 
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cal Bloch equations |l2j , and more recently Kimura who 
extended earlier work by Gorini to the strong coupling 
regime ,19j. Recently, there has also been considerable re- 
search activity on quantum Markov channels for two-level 
systems, motivated by their importance in quantum com- 
puting and communication. See, for instance, Ref. |2fij 
for a comprehensive analysis. A few simple higher di- 
mensional systems such as a three-level U-system with 
decay from two upper levels to a common ground state 
have also been studied, for instance by Lendi [T^ ]. 

In general, however, ensuring complete positivity of 
the evolution is often neglected for open systems with 
more than two (or degenerate) energy levels. For in- 
stance, the general expressions (6. A. 11) in [2l| for the 
relaxation rates ensure complete positivity only for two- 
level systems, as we shall show. For higher-dimensional 
systems additional constraints must be imposed if com- 
plete positivity is to be maintained. One reason for this 
neglect of positivity constraints is that, although the gen- 
eral form of the admissible generators for quantum dy- 
namical semi-groups is known, it can be difficult to verify 
whether a proposed dynamical law for an open system is 
consistent with positivity requirements. The main objec- 
tive of this paper is to address is issue. 

The paper is organized as follows. Starting with a 
purely phenomenological description of the interaction 
of an open system with its environment in terms of ob- 
servable population and phase relaxation rates — anal- 
ogous to the T\ and T-i relaxation times for a two-level 
system — we derive a general form for the dissipation 
superoperator in section [n] We then explicitly demon- 
strate with simple examples in section ITTT1 that the relax- 
ation rates cannot be chosen arbitrarily if the evolution 
of the system is to be physical in the sense that it sat- 
isfies complete positivity. In particular, we show that 
the phase relaxation rates for N > 2 are correlated even 
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in the absence of population relaxation, i.e., there exist 
constraints on the phase relaxation rates that are inde- 
pendent of population decay. To understand the nature 
of these constraints we express the empirically derived 
relaxation superoperator in Lindblad form fSection llVfl . 
and show that it can always be decomposed into two 
parts, one accounting for population relaxation and the 
other for pure phase relaxation processes ( Sect ion IVj l. 
This decomposition provides a general formula for the 
decoherence rates induced by population decay, which is 
consistent with physical expectations and evidence, and 
additional positivity constraints on the decoherence rates 
resulting from pure phase relaxation for N > 2. 

In section l\Tl we study the implications of these addi- 
tional constraints in depth for three and four-level sys- 
tems. In particular, we use the abstract positivity con- 
straints to derive explicit inequality constraints for the 
observable decoherence rates. Such explicit constraints 
are important from a theoretical and practical perspec- 
tive because they allow us to make concrete, empirically 
verifiable predictions about the decoherence rates and the 
dynamics of the system, as we show in section rVlll for sev- 
eral common, generic three and four-level systems such as 
ladder, A, V and tripod systems, and transitions between 
doubly degenerate energy levels. Experimental data con- 
sistent with positivity constraints would be significant 
and validate the chosen model for the open system dy- 
namics. On the other hand, if the observed relaxation 
rates for a system do not satisfy the constraints required 
to ensure complete positivity, it would be a strong indi- 
cation that the model used is not sufficient to properly 
describe the dynamics of the system. This does not nec- 
essarily mean that the model is useless; it might well be 
adequate for some purposes, but there will be cases where 
the model makes unphysical predictions and better mod- 
els, consistent with physical constraints, are required. 



II. QUANTUM LIOUVILLE EQUATION FOR 
DISSIPATIVE SYSTEMS 



density operator p based on the degrees of freedom of 
S, and describe its non-unitary evolution by amending 
the quantum Liouville equation to account for the non- 
Hamiltonian dynamics resulting from the interaction of 
S with the environment E. 

In this paper we restrict our attention to the (com- 
mon) case where the effect of the environment E leads 
to population and phase relaxation (decay and deco- 
herence, respectively) of the system S, and ultimately 
causes it to relax to an equilibrium state. To clearly de- 
fine what we mean by the terms population and phase 
relaxation, note that given an iV-dimensional quantum 
system we can choose a complete orthonormal basis 
{\n) : n — 1,2,..., N} for its Hilbert space and expand 
its density operator with respect to this basis: 
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(2) 

Although we can theoretically choose any Hilbert space 
basis, physically there is usually a preferred basis. Since 
the interaction with the environment usually causes the 
system to relax to an equilibrium state that is a statis- 
tical mixture of its energy eigenstates [25|, it is sensi- 
ble to choose a suitable basis of (energy) eigenstates of 
the system for modelling the relaxation process. In this 
setting the diagonal elements p nn in expansion Eq. (0) 
of p determine the populations of the (energy) eigen- 
states |ti), and the off-diagonal elements p nn > (n ^ n') 
are called coherences, since they distinguish coherent 
superpositions of energy eigenstates |W) = X^=i c nl n ) 
from statistical (incoherent) mixtures of energy eigen- 
states p = J2 n =i w n\n)(n\. 

Population relaxation occurs when the populations of 
the energy eigenstates change, typically due to sponta- 
neous emission or absorption of quanta of energy at ran- 
dom times. To account for population relaxation as a 
result of the interaction with an environment we must 
modify the system's quantum Liouville equation Q to: 



The state of an A^-level quantum system is usually rep- 
resented by a density operator p acting on a Hilbert space 
Tt. If the system is closed then its evolution is given by 
the quantum Liouville equation 



ih-p(t) = lH, P (t)l 



(1) 



where H is the Hamiltonian. Formally, the dynamics of 
an open system S that is part of a closed supersystem 
S + E (possibly the entire universe) is determined by 
the Hamiltonian dynamics (Q of S + E, and the state 
of the subsystem S can be obtained by taking the par- 
tial trace of the entire system's density operator ps+E 
over the degrees of freedom of the environment E. Of- 
ten however, the evolution of the (closed) super-system is 
unknown or too complicated and we are interested only 
in the dynamics of S. It is therefore useful to define a 
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P(*)])nn~y^ 7knPnn(t) + ^2 JnkPkk(t), 
k^n k=£n 

(3) 

where H represents the Hamiltonian dynamics of S, and 
7/t„ is the rate of population relaxation from state \n) to 
state \k), which depends on the lifetime of state \n), and 
in case of multiple decay pathways, the probability for the 
particular transition, etc. The jkn are thus by definition 
real and non-negative. Population relaxation necessarily 
induces phase relaxation, and we will later derive explicit 
expressions for the contribution of population relaxation 
to the phase relaxation rates. 

In general, phase relaxation occurs when the inter- 
action of the system with the environment destroys 
phase correlations between quantum states, and thus 
converts coherent superposition states into incoherent 
mixed states. Since coherence is determined by the off- 
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diagonal elements in our expansion of the density op- 
erator, this effect can be modelled as decay of the off- 
diagonal elements of p: 

p kn (t) = -~([ff,K*)J)ta - r knPkn (t), (4) 

where Tkn (for k ^= n) is the dephasing rate of the tran- 
sition |fc) <-> \n). 

Hence, population and phase relaxation change the 
evolution of the system and force us to rewrite its quan- 
tum Liouville equation as: 

Kt) = -±lH,p(t)]+L D \p{t)], (5) 

where Lr)[p{t)] is the dissipation (super-)operator deter- 
mined by the relaxation rates. It is convenient to note 
here that the N x N density matrix p{t) can be rewrit- 
ten as an N 2 (column) vector, which we denote as \p(t))), 
by stacking its columns. Since the commutator \H, p{t)\ 
and the dissipation (super-)operator Lj)[p(t)] are linear 



operators on the set of density matrices, we can write |J5J 
in matrix form: 

J t lP(t)))= (-~L ff + L D )|p(i)» (6) 

where Lh and Td are N 2 x N 2 matrices representing the 
Hamiltonian and dissipative part of the dynamics, re- 
spectively. Comparison with equations (0 and (0J shows 
that the non-zero elements of To are: 

(^)(m,n),(m,n) = -^mn m ^ n 

(I j u)(m,m),(m',m') = +7mm m ^ ml (7) 

^D){m.m),(m,m) = ~ 2j 7fcm 

where the index (m, n) should be interpreted as m + 
(n - 1)N. r mn = r nm implies {Ln)( m ,n),(m,n) = 

For a three-level system subject to population and 
phase relaxation, for instance, equation J7J gives a dissi- 
pation super-operator of the form: 
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(8) 



where jkn and Tkn are the population and phase relax- 
ation rates, respectively. 



III. PHYSICAL CONSTRAINTS ON THE 
DYNAMICAL EVOLUTION 

Although J7J gives a general form for the dissipation 
superoperator of a system subject to population and 
phase relaxation, not every superoperator of this form 
is acceptable on physical grounds since the density op- 
erator p(t) of the system must remain Hermitian with 
non-negative eigenvalues for all t > 0, and its trace must 
be conserved |26| . It is easy to see that the relaxation 
parameters in Q cannot be chosen arbitrarily if we are 
to obtain a valid density operator. For instance, it is well 
known in quantum optics that a two-level atom with de- 
cay |2) — > |1) at the rate 712 > also experiences dephas- 
ing at a rate Ti2 > ^712 since the coherence P12 must de- 
cay with the population of the upper level in order for p(t) 
to remain positive, consistent with the constraints on the 



relaxation rates for two-level systems derived in |lg, [19| . 

In higher dimensions we also expect population relax- 
ation from state \n) to |fc) at the rate 7fc„ to induce de- 
phasing of this transition at the rate Tkn > \lkn- How- 
ever, for N > 2 the situation is more complicated. First, 
a single random decay \n) — > \k) due to spontaneous 
emission, for instance, may affect other transitions in- 
volving the states or |n). This is perhaps not too 
surprising but since it is a crucial motivation for the fol- 
lowing sections, we shall consider two concrete examples. 

First, consider a three-level system subject to decay 
\2) — > |1) at the rate 712 but no other relaxation. Sup- 
pose, for instance, we follow formula (6. A. 11) in [2l| and 
set Ti2 = 7;7i2 and take all other relaxation rates to be 
zero. Then, assuming H = for convenience, the solu- 
tion of Eq. (JSJ for this dissipation super-operator leads 
to the density matrix 

//7ii + (l-e-*^ 2 )p 2 2 e-^^pxi p 13 \ 
P(t) = e-*W2 p21 e-^paa P23 , 

\ P31 P32 P33 / 

(9) 
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which in general is not positive for t > 0. For example, 
the superposition state 

p(0) = |*)(*|, |*) = -^(1,1, If (10) 

evolves under the action of this dynamical generator to 
a "state" p(t) which has a negative eigenvalue (i.e., neg- 
ative populations) for alH > as shown in Fig. ^ and is 
thus physically unacceptable. 

Furthermore, population relaxation is not the only 
source of constraints on the decoherence rates for N > 2. 
A perhaps more surprising observation is that even if 
there is no population relaxation at all, i.e., 7fc„ = for 
all k, n, and the system experiences only pure dephasing, 
we cannot choose the decoherence rates T^n arbitrarily. 
For example, setting T 12 ^ and T 2 3 = T 13 = for our 
three-level system gives 

/ Pn e~ ri2t pi2 P13 \ 

p(t) = e" V ^ t p21 P22 P23 ■ (11) 

V P31 P32 P33 J 

Choosing p(0) as in Eq. we again obtain a density 
operator pit) with negative eigenvalues (see Fig.^}. This 
shows that there must be additional constraints on the 
decoherence rates to ensure that the state of the system 
remains physical. 



IV. STANDARD FORM OF DISSIPATION 
SUPER-OPERATORS 

Significant progress toward solving the problem of find- 
ing dynamical generators for open systems that ensure 
complete positivity of the evolution operator, and hence 
positivity of the system's density matrix, was made by 
Gorini, Kossakowski and Sudarshan [16j who showed that 
the generator of a quantum dynamical semi-group can be 
expressed in standard form 

L[p(t)\ = -i[H,p(t)] 

N 2 -l 

+ 2 E akk'(lV k p(t),v£,J + lV k ,p(t)v£i) 

k,k' = l 

(12) 

where H is the generator for the Hamiltonian part of the 
evolution and the Vk, k = 1,2,..., TV 2 — 1, are trace- 
zero, orthonormal operators (Vk,Vk>) := Tr(V^Vfc/) = 
Skk' that together with V^ 2 — form a basis for 

the system's Liouville space. Furthermore, the resulting 
evolution operator is completely positive if and only if the 
coefficient matrix a = iflkk' ) is positive. 

Noting that a positive matrix {auk' ) has real, non- 
negative eigenvalues jf. and can be diagonalized by a 
unitary transformation, we obtain the second standard 
representation of dissipative dynamical generator, which 



was first derived (independently) by Lindblad |14| : 
L[p{t)\ = -i[£T,p(t)] 

+2 E 7*(n*p(t),4i + u k ,p(t)Aii). 
fc=i 

(13) 

Yet, although the general expressions Ijl2(l and ljl3|l have 
been known for more than two decades, it is often un- 
known whether a proposed generator for the dissipa- 
tive dynamics for a particular model is completely posi- 
tive, and some common dissipative generators have been 
shown not to satisfy this condition, as in the case of the 
Agarwal/Redfield equations mentioned earlier. In part 
this may be due to the fact that it is often very difficult 
in practice to express phcnomenologically derived dissi- 
pation generators in either of the two standard forms, and 
hence to verify if a proposed generator satisfies complete 
positivity. 

However, given a matrix representation for the relax- 
ation super-operator of the form J7J , which was derived 
from a purely phenomenological model based on observ- 
able decay and decoherence rates, we can express it in 
standard form (|12f> and transform abstract positivity re- 
quirements into concrete, easily verifiable constraints on 
the empirically observable relaxation rates. For this pur- 
pose we need a basis {Vk} for the Liouville space of the 
system. A canonical choice is to define N — 1 diagonal 
matrices 

V(m,m) = , ^ = ( E Gss ~ 171 e ™+l,m+l ) (14) 

for m = 1, 2, . . . , N — 1, as well as N 2 — N off-diagonal 
matrices 

V(m,n) = e mn, m ^ 11, 171, U = 1, 2, . . . , N. (15) 

where e mn is an N x N matrix whose entries are zero 
except for a one in mth row, nth column position. It 
is quite easy to verify that the N 2 — 1 operators V( m>n ) 
thus defined are trace-zero N x N matrices that satisfy 
the orthonormality condition Tr^VkVu,) — Skk' for any 
dimension N. 

Having defined the basis operators Vk, we can now 
compute the generators 

Lkk'W)] = \ (iV kP {t),V k \\ + iV k ,p{t)V k \\) (16) 

of the dissipation super-operator (|12fl with respect to 
this basis, where k, k' = 1,2,..., iV 2 — 1. Recalling that 
L kk '[p(t)] is equivalent to Lkk'\p(t)}}, where each hkk' is 
an N 2 x iV 2 matrix and \p{t))) is an iV 2 -column vector, 
we note that any (trace-preserving) dissipation superop- 
erator ho can be written as a linear combination of these 
dissipation generators: 

N 2 -l 

Ll> = E a kk'^kk'- (17) 
k,k' = l 
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FIG. 1: Eigenvalues of Eq. ©, left, and Eq. JTTJ, right, for p(0) as in Eq. (O 



To compute the coefficient matrix a = (akk 1 ) we can 
rewrite the TV 2 x N 2 matrices L^k' and L^> as column 
vectors Lkk' and of length iV 4 , and a as a column 
vector a of length (iV 2 — l) 2 , and solve the linear equa- 
tion Ld — Aa where A is an iV 4 x (iV 2 — l) 2 matrix 
whose columns are the hkk' ■ This matrix equation has 
a solution for any trace-zero Liouville operator since the 
columns of A span the space of trace-zero Liouville op- 
erators of the system. This procedure allows us in prin- 
ciple to express any (trace-preserving) dissipation super- 
operator in standard form (|12J) . and to verify whether 
it generates a completely positive evolution operator by 
checking the eigenvalues of the coefficient matrix (akk' ) • 
However, in practice this is not very efficient, especially 
for large N. Instead, we would like to be able to express 
the coefficients akk' directly in terms of observable relax- 
ation rates. This is the aim of the following sections. 

V. DECOMPOSITION OF RELAXATION 
SUPEROPERATOR 

We now use Eq. I|17|) to show that the relaxation super- 
operator of any iV-level system subject to both popula- 
tion and phase relaxation processes can be decomposed 
into a part associated with population relaxation pro- 
cesses and another accounting for pure decoherence. To 
this end, we introduce two types of decoherence rates: 
and T d nn for decoherence due to population re- 
laxation and pure phase relaxation (dephasing), respec- 
tively, and require that T mn — T P nn + Tf ml . 

If we have population relaxation \n) — > \m) at the rate 
Tmn > for m,n = 1,2, ...,N (with j mm = 0), then 
setting akk — Iran for k = m + (n — l)N and akk' = 
otherwise in lfT7|) leads to a dissipation superoperator 

N 

= ^ ' 7mn^(m,n),(m,n) ■ (1$) 
m,n=l 



Inserting Eq. i|16fl for h ( m ,n ),(m,n) with k = m+(n — 1)N 
and Vk as in Eqs (|14|l -l|15 |) . we obtain 

(I J £))("i,m),(m,m) = ~ Etl,fc^m Ttmi 
(^ J I))( m 5 m ) 1 ( m/ ' m/ ) 7mm' i TTi ^ Tfl 

(L£,)( m ,„),(m,„) = ~\ 2fc=l(7fem + 7fen), m^n 

which agrees with the general form Eq. (0) of the relax- 
ation super-operator, yields the correct population relax- 
ation rates, and suggests that the dephasing rates due to 
population relaxation are given by: 

1 N 

fe=l 

i.e., that the decay-induced decoherence of the transition 
between states \n) and \m) is one half of the sum over 
all decay rates from either of the two states \m) or \n) to 
any other state. Finally, inserting 

1 N 

r mn = T d mn + - 2j(7fe m + 7fcn), m^n (20) 
fc=i 

into Eq. 10 and solving Eq. (|17|1 shows that the dis- 
sipation super-operator hp of the system decomposes, 
L D = L P D + hp, with L P D given by Eq. (JTSJ and 

N-l 

L^) ^ ® (m ,m) ,(m' ,m J (m ,m) ,(m' ,m' ') ■ (^-0 

m,m' = l 

Thus, given L^> and the population relaxation rates j m n 
of the system, we can compute \j d D = hp — \P D and de- 
termine the coefficients b mm > := a( m . m ),(m'.m') m JUt by 
rewriting the super-operators \j d D and h( m ,m) ,(m',m') as 
column vectors I and Ik, respectively, defining a matrix 
B whose columns are given by Ik, and setting b = B~ l l 
where B~ x denotes the pseudo-inverse of B. However, 
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note that the matrix B only has (N — l) 2 instead of 
(iV 2 — l) 2 columns, and we can eliminate all zero rows. 
The resulting coefficient vector b can be rearranged into 
an (N — 1 ) x (N — 1) coefficient matrix b — (b mm >) that 
depends only on the pure dephasing rates Further- 
more, the requirement of positivity of the coefficient ma- 
trix (flkk' ) in Eq. (|f 2fl now reduces to the (much simpler) 
requirement that the (N — 1 ) x (N — 1 ) matrix (b mm i) 
be positive semi-definite. 

It is important to note that our formula l|19|) for the 
contribution of population relaxation to the overall deco- 
herence rates, obtained solely by imposing the physical 
constraint of complete positivity on the evolution of the 
system, a gree s with the expressions given, for instance, 
by Shore [22] for the general Bloch equations of iV-level 
atoms subject to various dissipative processes, but our 
general expression for the dissipation super-operator cov- 
ers systems subject to both population decay and pure 
dephasing processes, and implies the existence of non- 
trivial constraints on the pure dephasing rates of the sys- 
tem for N > 2. In the following sections, we shall analyze 
these constraints in detail for N = 3 and N = 4. 



VI. CONSTRAINTS ON THE PURE 
DEPHASING RATES 

A. Three-level Systems 

Expanding the relaxation super-operator JSJ) for a 
three-level system with respect to the basis 



we obtain 



V(2,2) 



= 



- e 2 2), 

+ e 22 - 2e 33 ) 

m, n — I, 2, 3, m =/= n 



where e mn is the 3x3 matrix whose entries are zero 
except for a I in the mth row, nth column position — 
which corresponds to the canonical basis (|14|) - (|15[1 for 

N = 3 — yields an 8 x 8 coefficient matrix (akk 1 ) whose 
non-zero entries are akk = Imn for fc = m + 3(n— 1) and 
m ^ n, as well as 0( m , m ),( m / )m /) = b mm >, where 



6n 
hi 

&22 



-Ix/3 



|V3 
4V3 



3 
4 



3 
4 
3 
4 
3 
4 



T d 

1 12 

T d 

1 23 

rf 3 

V r 23 



bn 



12 



&22 = (-rf 2 + 2rf 3 + 2r^ 3 )/3 (23) 
bi2 = b 21 = (r d 13 -r d 23 )/Vs 

Therefore, the coefficient matrix (akk' ) of the relaxation 
superoperator ho will be positive semi-definite if and 
only if 7„ m > and the real symmetric 2x2 matrix 



611 612 

&21 &22 



has non-negative eigenvalues. The second condition is 
equivalent to 



bn + 622 > 0, 611622 > b\ 2 . 
Substituting ill' Ml) these conditions become 



(24) 



< r? 2 + rf 3 



r^ 3 <2 



rf ? r d 



2 X 13 



1 12 1 23 



1 13 1 23 



(25) 

This double inequality provides upper and lower bounds 
for the pure dephasing rates of the system. As expected 
we can also show that each relaxation rate must be non- 
negative since the second inequality is of the form [a + 
(b + c)] 2 < 4 [a (6 + c) + be], which can be rewritten as 
a 2 + (b — c) 2 < 2a(b + c), and the first inequality implies 
a + b + c> 0, which shows that a > and b + c > 0, and 
by symmetry of a, 6, c requires that Tf 2 , and T d 3 be 
non-negative. Hence, we can rewrite (|25f) as follows: 



(26) 



( - vT c r < r a < (vt 6 



(22) where {a, b, c} is any permutation of {12, 13, 23}. 



Note that the choice T d 2 > and rf 3 = = 0, 



yd 

1 12 



which corresponds to the second example in section IIIII 
if there is no population relaxation, clearly violates 12611 , 
which explains why it results in non-physical evolution. 
We also see that allowed choices include, for instance, 
^ 3 > and rf 3 = 0. 
In general, H26fl shows that pure dephasing in a three- 
level system always affects more than one transition. Fur- 
thermore, if two of the pure dephasing rates are equal, say 
T d , then the third rate must be between and 4r d . For 
instance, consider a triply degenerate atomic energy level 
with basis states \m = 0) and |m = ±1). If rl 10 = Fq j 
then < h d _ 11 < ATq 1: i.e., the decoherence rate of the 
transition between the outer states can be at most 4Fq x . 
But note in particular that it could be zero even if the 
decoherence rate between adjacent states is non-zero. 



and the pure dephasing rates T d nn are defined by (|20|l . 
Noting that the pseudo-inverse of the matrix B is 



Four-level Systems 







/ 5 i 

^ L ^ L 

2v/3 \/3 2\/3 V3 

^_ ^ ^ x - 

2V3 2V3 V3 V3 

_1 1 1 11 1 / 

6 3 6 3 3 3 ' 



If we expand the relaxation superoperator h]j for a 
four-level system as discussed in section |3 with respect 
to the standard basis 114|) - (|15(l . we again obtain a coef- 
ficient matrix (ct( m ,n),(m',n')) whose non-zero entries are 
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a( m ,n),(m,,n) = Iran for TO ^ n, as well as a(ro,m),(m',m') = 

&mm' with bn, b 12 , b 21 and 6 22 as in (J23j, 6 31 = 6 13 , 
^32 = ^23 and 



as well as 



b 23 
b 3 3 



V6(-rf3 + 3rf 4 + rf 3 -3rf 4 )/i2 

V2(-2Tf 2 + Tf 3 + 3T d 14 + T d 23 + 3T d 24 - 6rf 4 )/12 

(-rf a - r? 3 + 3rf 4 - r^ 3 + sr^ + srf 4 )/6 

(27) 

Since the reduced coefficient matrix b = {b mm ' ) is a real, 
symmetric 3x3 matrix, necessary and sufficient condi- 
tions for it to be positive semidcfinite are |23| : 



fell > 0, 



bllb 22 > b\ 2 , 



det(6) > 0. 



(28) 



The first two of these conditions are equivalent to (|24H . 
Thus, the pure dephasing rates for a four-level sys- 
tem must satisfy (|25H and (|26|) . and the new constraint 
det(6) > 0, or: 



bub 22 b 33 



2b 12 b 13 b 23 > &11&23 



b 22 b 2 3 + b 33 b 2 2 . (29) 



Unfortunately, inserting (|23|l and l)27|) into this inequality 
does not yield a nice form for the additional constraint. 

We can obtain a more symmetric form of the con- 
straints by choosing a slightly different operator basis: 



V (2,2) 

V {3,3) 

V! , 
(m,n) 



i(en - e 22 + e 33 - e 44 ) 
|(e u - e 22 - e 33 + e 44 ) 
2 ( e n + e 22 - e 33 - e 44 ) 
m, n = 1, 2, 3, 4, 



(30) 



2 



TO 7^ 71 



The V7 m „•) are trace-zero matrices that differ from the 
standard operator basis only in the choice of the diagonal 
generators and also form an orthonormal basis for the 
trace-zero Liouville operators of the system. However, 
expanding L d with respect to this basis (|30|l gives a more 
symmetric coefficient matrix b' with non-zero entries: 



(31) 



where T? ot = \ £* =2 YT m =i T tn is half the sum of all 
pure dephasing rates. Since the eigenvalues of the coeffi- 
cient matrix are independent of the operator basis, b' has 
the same eigenvalues as b. Furthermore, necessary con- 
ditions for b' to have non- negative eigenvalues are |27|: 





b'u 


pd 
— 1 tot 






b'22 


- Y d 

— 1 tot 


-(rf 4 + r^ 3 ) 




bk 


- r d 

— 1 tot 


-(r? 2 + i1 4 ) 


b' 12 


= b' 21 


= (rf 2 


-rf 4 )/2 


b'u 


= b 3 i 


= (rf 4 


- r1 3 )/2 


b' 2 . 


= b'32 


- (rf 3 


-rl)/2 



hi > 0, 6 22 > 0, 
bnb 2 2 > b\ 2l 6n6 33 > b [ 



2 

13; 



633 > 0, 

^22^33 > b 



23- 



Inserting (JHU) into (j23) and ||23J) yields 



1 13 
1 14 

1 12 



24 



< rf 2 



^23 — ^12 

pd < r rf 

1 34 ^ 1 13 



r? 4 
rf 4 



^24 
1^23 



L 34 

T d 
1 34 

r d 

J- 9/1 



(32) 
(33) 



(34) 



(r? 4 
(rf 2 
(rf 2 



1 23 1 13 



r 3 d 4 



1 13 
pd 
L 14 



1 24J 

- r d } 2 
1 24) 



r d 2 3 ) 2 



4 pd pd 
4i 12 1 34 
^pd pd 
41 14 1 23 
^ipd pd 
* L 13 1 24 



(35) 



which can be written simply as 

\b-c\ < a < b + c, (b- 



c) 2 < Axy 



(36) 



if {a,b,c} is a permutation of the set {rf 2 + r 34 ,Tj 



pd pd 
1 24; 1 14" 



~-T 23 } and we let x and y be the summands of a, 



e.g., if a = rf 2 Tx 34 
a = x + y shows in particular that < x + y and < Axy 
and thus x, y > 0, from which we can conclude especially 
that T d nn > for all to, n. 

In certain cases these constraints can be simplified. For 
instance, if the pure dephasing rates for transitions be- 
tween adjacent states are equal, i.e., Tf 2 — T 23 = T 34 = 
rf , as one might expect, for example, for a system con- 
sisting of the basis states of a four-fold degenerate energy 
level, and we set r 2 = |(rf 3 + T^ 4 ) then (gH) yields the 
following bounds on the decoherence rate Tf 4 : 



T^ 4 then x = Tf 2 and y = Tf 4 . Sett 



mg 



max{2r 



3rf,rf 



2r^, 0} < T u <T d + 2T^ (37) 



and combined with the second inequality of (|35|l we ob- 
tain <Ti < 4rf, and thus Tf 4 < 9Tf and rf 3 ,rf 4 < 



8Tf. If we further assume that the pure dephasing rates 
for transitions between next-to-nearest neighbor states 
are equal as well, i.e., Tf 3 = Tf 4 — T 2 , then we obtain 
rf 3 = T 24 < ATf. Hence, for a system whose pure de- 
phasing rates depend only on the "distance" between the 
states, the former are bounded above by 

k) 2 T' 



L nk <(n~ nj x l5 

where Tf = T d ln+1 is the dephasing rate for transitions 
between adjacent sites. 

In this special case we can compare the constraints ob- 
tained from our necessary conditions with the necessary 
and sufficient conditions (|2*5|l . Inserting T d k — rtu, 
into (|3*TJl yields 



(38) 



b'u - 

b'22 = 

b'33 = 

b[ 2 = b^ = 



|(3rf-2r^ + r^) 
|(rf + 2r^-rf) 
(-rf + 2r 2 i - 3rf ) 



21 



b'=b' 



J 13 

b'23 = b' 3 



2 


2^3 




(39) 



T d ) 



and the necessary and sufficient conditions l|28|) become: 



r 3 > 2r^ - 

rf < 2rf 
rfrf > (rf 



3rf 
-rf 

-T d ) 2 



(40) 



If Tf = then T d = 2T^ = 0. Otherwise, we can multiply 
the second inequality by Tf and combine it with the third, 
which leads to rf (2^ + Tf ) > (rf - r 2 ) 2 and simplifies 
to T 2 < 4rf . Inserting this result in the first inequality 
gives T 3 < 9rf, i.e., the necessary conditions (|38|) are 
also sufficient. 
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VII. EXAMPLES AND DISCUSSION 2 



We now apply the results of the previous sections to 
several types of generic three and four-level atoms. The 
objective in each case is to derive a proper relaxation su- 
peroperator, which is consistent with both experimental 
data and positivity constraints, and to discuss the impli- 
cations of the latter constraints. Though the emphasis is 
on atomic systems, the results generally apply to molec- 
ular or solid state systems with similar level structures 
as well. 




2 



A. Generic three-level atoms 

Let us first consider the general case of a generic three- 
level system subject to arbitrary population and phase 
relaxation processes. Let 7 m „ denote the observed rate 
of population relaxation from state \n) to state \m) for 
m,n — 1,2,3 and m ^ n, and let r 12 , r 23 and T 13 be 
the observed decoherence rates for the 1 «-> 2, 2^3 and 
1 <-> 3 transitions, respectively. Then the pure dephasing 
rates of the system according to Q20fl are: 



1 12 

1 13 

-pcZ 
1 oq 



Ti2 — (721 
Ti3 — (721 
r 2 3 — (712 



731 
731 
' 732 



- 712 
713 
713 



732)/2 

72 3 )/2 
72 3 )/2 



(41) 



These dephasing rates must be non-negative and satisfy 
the inequality constraints (|26|l for the evolution of the 
system to be completely positive. Experimental data for 
the observed relaxation rates that fails to satisfy these 
conditions should be considered a reason for concern, and 
might suggest that the physical system under investiga- 
tion cannot be adequately modelled as a three-level sys- 
tem, for instance. 

If the dephasing rates do satisfy the necessary con- 
straints then a physically valid representation of the re- 
laxation superoperator L d for the system in terms of the 
observed relaxation rates is 

LdW)] = ]T lm nLP mn [p(t)]+SiL d Mt)}+8 2 L d Mt)}, 

(42) 

where the elementary relaxation terms are 

VtnJptf)] = lV {m , n) p(t),V^ tn) } + lV {m!n) ,p(t)Vl mtn) j 

2L d m [p{t)] = lA m p{t),A m } + lA m ,p(t)A m j 

(43) 

with V( m .„) as in and the diagonal "pure dephasing" 
generators 



Ax 



A 2 = 



1 



y / 2x(x - Ai) 
-1 

^2x(x + A 1 ) 



n/3A 2 F ( i,i) + (x - Ai)F (2 , 2) 
\/3A 2 V (M) - (x + AijVfa.a 



FIG. 2: Three-state atoms: ladder system (left), A-system 
(top right) and V-system (bottom right) with arrows indicat- 
ing population decay pathways. 



and the "effective dephasing" rates are 

6 1/2 = (Ti 2 +Ti 3 +Tl 3 ±x)/3, 



(45) 



where x = ^/Af + 3A|, Ai = 2hi 2 - Ti 3 - r 23 and 

a 2 = — r 23 . 



1. Ladder configurations 

For a three-level atom in a ladder configuration where 
the main source of population relaxation is spontaneous 
emission from the excited states to a stable ground state, 
as shown in Fig. [21 we simply set 721 =731 = 732 = 
in l(4*T|) and l(4*2")l . respectively, to obtain the correct pure 
dephasing rates 



-pd 
1 12 


= r 12 


" 712/2 


j^d 
1 13 


= r 13 


- (713 + 723)/2 


-pd 
1 23 


= 1^23 


- (712 + 713 + 723)/2 



(44) 



(46) 



and the corresponding relaxation superoperator hjj . The 
decay-induced decoherence rates in this case satisfy 
the interesting equality r^ 2 + r^ 3 = r 23 . 

This is another way of seeing that, if 712 > as in 
example 1 considered in section ITTT1 then r 23 must be at 
least 712/2 — recall that we showed explicitly that the 
naive guess Ti 2 = 712/2 and Ti 3 = T 23 = leads to 
non-physical states with negative eigenvalues. 

r i3 = 0, on the other hand, is possible even if 712 > 
provided that state |3) is stable. It is interesting to note, 
however, that Ti 3 = always implies Ti 2 = T 23 . If there 
is no pure dephasing then this is obvious since r^ 2 = r 23 , 
but it is true even if there is pure dephasing, since Ti 3 = 
implies Tf 3 = and thus the inequality constraint l|26|) 
for the pure dephasing rates implies Tf 2 = r 23 . 
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2. A systems 

For a A system for which only the decay rates 712 and 
732 are non-zero as shown in Fig. |3 the pure dephasing 
rates are: 

r?2 = ri2-( 7 i2 + 73 2 )/2 

I*, = r 23 -(7i2 + 73 2 )/2 (47) 

r? 3 = r 13 . 

Moreover, if the lifetime of the excited state is 7 -1 
and the system is symmetric, i.e., 712 = 732 = 7/2 
and Ti2 = r 2 3, as is often the case, then we have 

r 12 = r 23 = r d + 7 /2. ifr d = othcnrf 3 = o due to (ESJ. 

Otherwise, setting Tf 3 = aT d gives A 1 = 2(1 - a)T d , 
A 2 = -(1 - a)T d , x = 2(1 - a)T d , and thus the relax- 
ation superoperator is simply 

L D [p(t)\ = h/2){LUp(t)]+LlMt)}} 

+ [{A-a)T d /i\L d [p{t)]+aT d L d [p{t)] 

where the diagonal pure dephasing generators 144(1 are 

A 1= A\ = (-V3V (lil) +V(2,a))/2 

A 2 =A\ = (V (1<1) + \/3V (2j2) )/2. (48) 

Positivity requires 7 > 0, 4 — a > and a > 0, or 
< Tf 3 < 4T d , in accordance with $3$ and previous 
observations. 

3. V systems 

Similarly, for a V system for which only the decay rates 
721 and 723 are non-zero as shown in Fig. |2 the pure 
dephasing rates are simply: 

r i2 = r i2 - 72l/2 

r& = r^ 3 - 723/2 (49) 
r? 3 = r1 3 -( 7 2i + 723)/2 

If the system is symmetric, i.e., both excited states have 
the same lifetime, 721 = 723 = 7, and Ti2 = T23 then 
setting rf 3 = aT d leads to the relaxation superoperator 

L D [p(t)\ = l{L p 2 Mt)]+L p 2 Mt)}} 

+ [(A-a)T d /3]L d 1 [p(t)}+aT d L d 2 [p(t)} 

with L mn [p(t)] and Li[p(t)] as defined in l|43l) and the 
generators A\, A 2 as in (@BJ). 

4- Comparison of A and V systems 

If the excited states have the same lifetime 7 _1 and the 
pure dephasing rate T d for transitions between the upper 
and lower states is the same for the A and V configura- 
tion, then we have Ti2 = T23 = 57 + T d in both cases, 



i.e., the overall decoherence rate for transitions between 
ground and excited states is the same for both configura- 
tions. The main difference, as expected, is the decoher- 
ence rate of the 1^3 transition, which is Ti3 = Tf 3 for 
the A-system, and r i3 = 7 + Tf 3 for the V^-system. 

Thus, if pure dephasing of the 1 <-> 3 transition is negli- 
gible then it will remain decoherence free for the A system 
but not for the V system. However, if pure dephasing is 
taken into account then the transition between the de- 
generate ground states of the A system may not be de- 
coherence free, and comparison of the decoherence rates 
for both systems, T^ 3 = (rf 3 ) A and r% = 7 + (rf 3 ) y , 
shows that r^ 3 could theoretically even be greater than 
r]^ if the pure dephasing rate of the transition between 
the degenerate ground states was greater than the decay 
rate 7 plus the pure dephasing rate Tf 3 for the V system. 



B. Generic four-level atoms 

Again, we will first consider the general case of a 
generic four-level system subject to arbitrary population 
and phase relaxation processes. Let 7 mn denote the ob- 
served rate of population relaxation from state \n) to 
state \m) for m,n = 1, 2, 3, 4 and m ^ n, and T mn be the 
observed decoherence rates for the m <-> n, transitions, 
as usual. Then the pure dephasing rates of the system 
according to lj2"U|) are: 



1 12 


= r 12 - 


(721 " 


1- 731 - 


h 741 - 


1-712 - 


1- 732 " 


h 74 2 )/2 


j^d 
1 13 


= r 13 - 


(721 " 


1- 731 - 


1- 741 - 


1-713 - 


1- 723 - 


h 74 3 )/2 


j^d 
1 14 


= r u - 


(721 " 


1- 731 - 


1- 741 - 


I"7l4 - 


1- 724 - 


h 734)/2 


j^d 
1 23 


= T23 — 


(712 " 


1- 732 " 


1- 742 - 


1-713 - 


h 723 " 


H74 3 )/2 


j^d 
1 24 


— T 2 4 — 


(712 " 


h 732 - 


1- 742 - 


1-714 - 


h 724 - 


h 734)/2 


1 34 


= T23 — 


(713 " 


1- 723 " 


1- 743 - 


1-714 - 


h 724 " 


h 734)/2 



(50) 



These dephasing rates must satisfy the necessary and suf- 
ficient conditions (|28|l for complete positivity, and in par- 
ticular the inequality constraints (|34|) and (|35|l . If the 
data for a given system does not appear to satisfy these 
conditions then (unless the data is unreliable) it should 
be assumed that the system cannot be properly modelled 
as a four-level system subject to population and phase 
relaxation. As mentioned in the introduction, such mod- 
els are quite common and may still be adequate for some 
purposes but can lead to non-physical results such as 
states with negative eigenvalues, etc. 

If the dephasing rates do satisfy the necessary con- 
straints then a physically valid representation of the re- 
laxation superoperator Ld for the system in terms of the 
observed relaxation rates is 

3 

LdW)] = ImnULnW)] + E b ™L d mn {p(t)} 

m^n m,n=l 

(51) 
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where the elementary relaxation terms are 



2LZ 



lV(m,n)P\"J, r ( m ,n)U 

W)] = lv( mim)P (t),v( n>n) \ 



Wim.n)-, HV) " ( m ,n) 

(52) 



with V ( m .n) as defined in V (m ,„) = V( m n) for m ^ n, 
and coefficients b mn as in (|31[1 . Note that — unlike in 
the three-level case — we chose not to diagonalize the de- 
phasing superoperator since the general expressions are 
quite complicated and do not confer a significant compu- 
tational advantage. 



Transition between doubly-degenerate levels 




FIG. 3: Four-state atoms: transition between doubly degen- 
erate energy levels (top), tripod system (bottom left) and 
inverted tripod (bottom right) with arrows indicating popu- 
lation decay pathways. 



The results of the last paragraph apply, for instance, 
to a system consisting of two doubly degenerate energy 
levels subject to population relaxation as shown in Fig. El 
and general phase relaxation. The decay-induced deco- 
herence rates according to (fT^|) are: 



pp — pP 

1 12 — 1 23 

pp — pP 

1 14 — 1 34 

1 24 
F P 
1 13 



(712 + 732)/2 

(714 +734) /2 

(712 + 732 + 714 + 734)/2 

0. 



(53) 



Thus, if all decoherence is the result of population relax- 
ation processes then the transition between the ground 
states remains decoherence free. However, if there is 
pure dephasing then i|3rj|) implies (b — c) 2 < Axy for 



y with x = Tf 3 , y = P 



24 J 



r d 
1 12 



pd 
1 34 



a = 

c = Tf 4 +~rf 3 . Thus, x = ff 3 = is possible if and 
only if b = c. Since T p 4 + T 23 = T\ 2 + T 34 according to 
(|53[1 . this is equivalent to Ti2 -I- r 34 = Fi4 + 1^23. Con- 
versely, if Fi2 + r 3 4 7^ Ti4 + T23 then b c, and we have 
0<\b — c\<x + y and < (b — c) 2 < 4xy, which implies 
x > 0, i.e., Tf 3 > 0. 

Now suppose both excited states have the same Ti- 
relaxation time, i.e., the same spontaneous emission rate 
7, and the relative probabilities for the possible decay 
pathways are given by the absolute value of the Clebsch- 
Gordan coefficients of the transition. Then 712 = 734 = 
7/3 and 732 = 714 = 27/3, and the decay-induced deco- 
herence rates are Y\ 3 = 0, Y\ 2 — Y v 23 = = T 34 = 7/2 
and T 24 = 7, as one would reasonably expect. 

Furthermore, if the dephasing rates satisfy Tf 2 



i : 



rf and rf 4 = r d 23 



r d 

1 2> 



as one might expect 
due to symmetry for a typical system, then i|36|l im- 
plies especially (rf - T 2 ) 2 < 4rf 3 r 2 ! 4 . Thus, if we have 
rf 3 = then we must also have rf = T 2 , and conversely, 
if rf 7^ T 2 then rf 3 > 0, i.e., the transition between 
the two ground states can remain decoherence free only 
if rf = T^. This observation may seem trivial but it 
might be a convenient way of ascertaining if the transi- 
tion between the ground states is decoherence free or not 
by simply measuring the decoherence of the transitions 



between the ground and excited states, and the decay 
rates of the excited states. 

Moreover, if Ti 3 = 0, then we must have Tf 2 = rf 4 = 
= : T d according to our previous observations. 



1 23 



— 1 34 



If T d = as well then T^ 4 = due to ^ and we have 
Lf)[p(t)] = 0, i.e., no dephasing takes place. Otherwise, 
setting Y 2i = aT d leads to the simplified relaxation su- 
peroperator 

L D [p(t)} = (7/3){L{ 2 [p(t)}+L p 3 Mt)}} 

+(2 7 /3){i§ a K*)]+i? 4 K*)]} 

+ (4 - a)T d (lA lP (t), Axj + lA^pffiA!]) /4 

+aT d (lA 2 p(t), Ai\ + lA 2 ,p(t)A 2 j) /2 

(54) 



)/v / 2 and 
^ as defined in 13011 . 



^(2.2 



V (3,3)J 



< r° 4 < 4r d 



with A x = V[ xl y A 2 = (- 
LfnnW)} as defined in g3J} and V! 
Again positivity requires < a < 4, i.e. 
consistent with our previous observations, and thus pro- 
vides an upper bound of 4r d + 7 on the total decoherence 
of the transition between the upper levels. 



2. Tripod and Inverted Tripod Systems 

Another common type of four-level system is a tri- 
pod system, i.e., a transition between a triply degener- 
ate ground state and a non-degenerate excited state |4). 
With population relaxation due to spontaneous emission 
as indicated in Fig. [3J the decay-induced decoherence 
rates according to ifTHjl are: 



r? 2 
r? 4 



pp _ pp 

1 13 — 1 23 
pP _ pP 
1 24 — 1 34 





(714 +724 +734)/2 



(55) 



Assuming that the lifetime of the excited states is 7 1 
and all decay pathways are equally probable, we obtain 
7i4 = 724 = 734 = 7/3 and T^ 4 = T P 2A = Tf 4 = 7/2, as 
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well as 

l p d W)\ = hmmMt)} + ^Mt)] + HMt)}) (56) 

with as defined in (g3J). 

For comparison, the decay-induced decoherence rates 
(|19fl for an inverted tripod, i.e., a transition between a 
non-degenerate ground state |4) and a three- fold degener- 
ate excited state with population relaxation as indicated 
in Fig. are: 



1 12 — 


|(741 - 
2(741 " 


1- 742) 


pP _ 
1 13 — 


^ 743) 


pP _ 
1 23 — 


2(742 - 


1- 743) 


pP _ 
1 14 — 


2741 
9742 
2 743- 




pP _ 
1 24 — 




pP 

1 34 — 





(57) 



and assuming the lifetime of the excited states is 7 1 we 
obtain thus 741 = 742 = 743 = 7 and r? 4 = = T| 4 = 
7/2, T p 12 = rf 3 = r? 3 = 7 as well as 

L>(*)] = 7(LSxKt)] + LSa^t)] + L^t)]) (58) 

with L$n n [p(t)] as defined in 

Interestingly, the decay-induced decoherence rates for 
transitions between the upper and lower sublevels are 
the same for both systems. In absence of pure dephas- 
ing, the only difference between the two cases is that the 
degenerate subspace remains decoherence- free for the tri- 
pod system, while the decoherence rates for the inverted 
tripod are equal to the spontaneous emission rate 7 for 
the upper levels. This basic situation does not change 
very much if we add pure dephasing since the tripod and 
inverted tripod system are equivalent as far as pure de- 
phasing is concerned. However, there will be additional 
constraints, and we shall in particular study the case 

r? 4 = r^ 4 = rf 4 = r d and rf 2 = r^ 3 = r d , which 

one might expect to occur in many systems for reasons 
of symmetry. 

If T d = then the necessary conditions (|35[) can only 
be satisfied for T 2 = Tf 3 = 0, i.e., if there is no pure 
dephasing for transitions between the upper and lower 
sublevels then all pure dephasing rates are zero and all 
decoherence in the system must be due to population re- 
laxation. Hence, Lj-)[p(t)] — 0. Otherwise, set T 2 = aT d 
and Tf 3 = (3T d . Inserting these values into the coefficient 
matrix b' [Eq. 1)31(1 ] allows us to directly derive necessary 
and sufficient conditions for positivity of this matrix by 
computing its eigenvalues. 

det(6'/r d — XI) = 0- A) (A 2 - P X + q) 

where p = 3 + 2a — (3 and q = (4a — a 2 — (3)/2, shows 
immediately that the eigenvalues of b' are Ai = Tf 3 , and 
•^2/3 = (pi Vp 2 4q)/2. Hence, necessary and sufficient 
conditions for positive semi-definiteness of b' are Tf 3 > 0, 
p > and q > 0, and the dephasing superoperator can 



be written as: 

L d D [ P (t)} - Ti 3 (lA 1 p(t),A 1 l + lA 1 ,p(t)A 1 })/2 
+X 2 (lA 2 p(t),A 2 j + lA 2 ,p(t)A 2 })/2 
+X 3 ({A 3 p(t), A 3 j + fA 3 , p(t)A 3 j) 12 (59) 

where we have 

Ai = (V( 2<2) -V( 3<3) )/V2 

A 2 = [AaV( 1A) + + x)V( 2>2) -0 + x)V^ 3) ]/{2x) 
A3 = [^V( 1A) +0-x)V( 2t2) -0-x)V( 3t3) }/(2x) 

with V! m n ^ as defined in (|30|l . and a — a — I, $ = (3 — 1 

and x — y/8(a — l) 2 + — l) 2 . Furthermore, p > 
implies (3 < 3 + 2a and q > is equivalent to a(4 — a) > 
[3 > and hence implies < a, (3 < 4. 

VIII. CONCLUSIONS 

Starting with very basic assumptions we defined a sim- 
ple yet general relaxation superoperator, which should be 
adequate to describe a wide variety of open systems not 
too strongly coupled to their environment, solely in terms 
of experimentally observable quantities such as the pop- 
ulation relaxation and decoherence rates of the system, 
without imposing any restrictions on the types of popu- 
lation and phase relaxation that can occur. 

The advantage of a relaxation superoperator thus de- 
fined is that it can describe the observed dissipative dy- 
namics of the system in principle as accurately as we can 
measure the relaxation rates. Unfortunately, however, 
there are several problems with this approach. One is 
that it can lead to relaxation superoperators that do not 
preserve complete or even simple positivity, as we have 
explicitly shown for several examples. Since any violation 
of positivity effectively means negative or even non-real 
probabilities, this is serious problem. 

To avoid such problems one must impose constraints 
on the relaxation rates. We have analyzed the nature 
of these basic constraints by expressing our relaxation 
superoperator in the standard form for dissipative gen- 
erators of quantum dynamical semi-groups derived by 
Gorini, Kossakowski and Sudarshan. We have also shown 
that it is possible to decompose our generic relaxation su- 
peroperator into two distinct parts associated with pop- 
ulation relaxation and pure dephasing processes, respec- 
tively, and that the coefficients of the Kossakoswki gener- 
ators for the population relaxation part can be identified 
(usually uniquely) with the observed population relax- 
ation rates, the only restriction being the obvious one 
that the decay rates be non-negative. Most importantly, 
the expressions we obtain for the decoherence rates in- 
duced by population relaxation agree with similar expres- 
sions found in the literature. 

However, population relaxation is usually not the only 
source of decoherence. To account for other sources of 
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decoherence, we have introduced pure dephasing rates 
for each transition by subtracting the decoherence in- 
duced by population relaxation processes from the ob- 
served overall decoherence rates. These pure dephasing 
rates define the pure-phase-relaxation superoperator, and 
we express the coefficients of the Kossakoswki genera- 
tors for this part of the relaxation superoperator explic- 
itly in terms of these pure dephasing rates for three and 
four-level systems. These expressions, unlike the gen- 
eral expressions for the coefficients of the population re- 
laxation superoperator, are more complicated, and the 
requirement of complete positivity results in nontrivial 
constraints on the dephasing rates, which we have an- 
alyzed specifically for three- and four-level systems, al- 
though the same type of analysis can be performed for 
systems of higher dimension. 

Finally, we have applied these general results to study 
their concrete implications for several simple but com- 
monly used three- and four-level model systems such as 
A and V systems, tripod and inverted tripod systems and 
transitions between doubly degenerate energy levels. In 
each case we have attempted to make concrete predic- 
tions about inequality constraints and correlations of the 
decoherence rates demanded by the requirement of com- 
plete positivity which are experimentally verifiable. Such 
experimental tests of the constraints could be useful in 
various ways. Confirmation of the correlations would vin- 
dicate the semi-group description of the dynamics. On 
the other hand, violation of the constraints required by 
complete positivity would suggest that our model of the 
system is not really adequate to capture its real dynamics 
although it may still be useful for certain purposes. 

IX. EPILOGUE: POSITIVE MATRICES 

There has been a great deal of mathematical work on 
the properties of positive matrices, and papers such as 
"Some inequalities for Positive Definite Symmetric Ma- 
trices" [Siam J. Appl. Math. 19, 679-681 (1970)] by 
F. T. Man, "The Space of Positive Definite Matrices and 
Gromov's invariant" [Trans. Am. Math. Soc. 274, 239 
(1982)] by R. P. Savage, or "Positive Definite Matrices 
and Catalan numbers" [Proc. Am. Math. Soc. 79, 177- 
181 (1980)] by F. T. Leighton and M. Newman appear 
relevant to our problem at first glance. However, despite 



the connection to this work suggested by their titles, they 
really address rather problems. 

F. T. Man, for example, studies the problem of com- 
paring positive definite symmetric matrices, in particular 
answering the question under what conditions P > Q, i.e. 
P — Q positive, implies P 2 > Q 2 for positive definite sym- 
metric matrices P and Q. Unfortunately, these results 
are not applicable to our problem. However, they may 
be relevant for issues such as the comparison of density 
matrices. 

R. P. Savage considers the space of n x n positive 
definite matrices X with det(X) = 1 under isometrics 
X — > AXA T where A e SL(n,R) and shows that it has 
a collection of simplices preserved by the isomctries and 
that the volume of the top-dimensional ones has a uni- 
form upper bound. One could perhaps say that the den- 
sity matrices p of interest to us are positive, and that the 
dynamical Lie group of the system provides isometries of 
a sort, but our density matrices are positive matrices of 
trace one, which actually rules out det p = 1, since pos- 
itivity requires the eigenvalues to be non-negative and 
Tr(p) = 1 requires that the sum of these eigenvalues be 
one, which implies det(p) < 1 unless n = 1 and p = 1. 

Similarly, Leighton and Newman show that the num- 
ber ofnxn integral, triple diagonal matrices that are 
unimodular, positive definite and whose sub and super 
diagonal elements are all one, is the Catalan number 
in) K n 1)' wrucn i s an interesting mathematical re- 
sult but not relevant to our problem since our density 
matrices, although positive definite, are not usually tri- 
diagonal, and even if they were, the elements on the sub 
and super diagonal (the coherences) would have to be 
less than one for normalized density matrices. 

We thank the referee for bringing our attention to the 
rich mathematical literature on the subject of positive 
matrices. 
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